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We study a family of deterministic models for highway traffic flow which generalize 
cellular automaton rule 184. This family is parametrized by the speed limit m and 
another parameter k that represents a "degree of aggressiveness" in driving, strictly 
related to the distance between two consecutive cars. We compare two driving strategies 
with identical maximum throughput: "conservative" driving with high speed limit and 
"aggressive" driving with low speed limit. Those two strategies are evaluated in terms of 
accident probability. We also discuss fundamental diagrams of generalized traffic rules 
and examine limitations of maximum achievable throughput. Possible modifications of 
the model are considered. 

1. Introduction 

Transport phenomena in complex systems, in particular models of highway traf- 
fic flow, attracted much attention in recent years. Much of the effort was concen- 
trated on discrete stochastic models of traffic flow, first proposed by Nagel and 
SchreckenbergjL and subsequently studied by many other authors using a variety of 
techniqueso'ffiQ. In what follows, we shall study a family of purely deterministic 
traffic models. The only randomness comes from the fact that the initial configura- 
tion of cars is chosen at random. This family of models represents various driving 
strategies, either chosen by drivers (distance between cars) or externally imposed 
(such as the speed limit). Our "artificial highway" consists of an array of L cells. 
Each cell is either occupied by a single car or empty. Cars can move only to the 
right, and we assume periodic boundary conditions. Time is discrete. At each time 
step, each driver moves his car according to some specified rule. The evolution is 
synchronous, that is, all cars move at the same time. In the simplest model of the 
family, car at site i can either move to site i + 1 this site is empty, or not move if 
site i + 1 is occupied. Thus, the state of a given cell i depends only on cells 
i and i + 1. This model is equivalent to cellular automatonEI rule f84 if the state 
of an occupied site is 1, whereas the state of an empty site is 0. Under this rule, 
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which rule table is 

000 -> 0, 001 -► 0, 010 -> 0, Oil -> 1, 
100 -> 1, 101 -> 1, 110 -> 0, 111 -> 1, 

the density of l's is conserved, meaning that the number of "cars" does not change 
with time. Moreover, each car can move at most one site to the right during one 
time step, so the "speed limit" for this model is rn = 1. For a lattice of length L, 
the average speed (sum of all speeds divided by the number of all cars) is, therefore, 
always less or equal to m = 1. 

Let us now assume that we start from a random configuration of density p, and 
that cars move according to rule 184. It has been recently proved! that for large L 
and t, the average speed v t at time t equals 

G(p,t) ifp<|, 
vt={ l-o (1) 



P 



-Q(p,t) otherwise, 



e(,, t , = i-Ml^E. ( 2 , 

VTrf 



where 



Since lim^oo 0(/O, t) = 1, the average speed in the long time limit = 1 when 
the car density p is less than 1/2 and Voc = (1 — p)/p otherwise. In statistical 
physics terminology, the system exhibits a second order kinetic phase transition, 
where p is the control parameter, and v the order parameter. The critical point 
is at exactly p = 1/2, and at the critical point vt approaches its stationary value 
as t -1 / 2 . Away from the critical point, the approach is exponential, and it slows 
down asp comes closer to 1/2. In fact, when the lattice is finite (as in the real life), 
the behavior of this model is not significantly different from the L — oo case. It is 
sufficient to perform L/2 iterations in order to reach the stationary statei. These 
considerations are illustrated in Figure |], representing the average car speed versus 
car density after 500 iterations of a 1000-site lattice. Numerical simulations are 
compared to the theoretical expression of given by 

1 if /? < |, 

— otherwise. ^ 

P 

It is straightforward to generalize rule 184 to higher velocities. Let us consider 
two cars A and B such that A follows B, and denote by g the gap between them. 
Then, at the next time step, car A moves g sites to the right if <7 < m, and m 
sites otherwise. Such a rule will be referred to as lZ m ,i, where m denotes the speed 
limit. With this notation, TZi^, represents rule 184. The reason for the additional 
subscript "1" will become clear later. In this case, the phase transition occurs at 
p c = 1/(1 + to). In the free moving phase, i.e., below the critical point, all cars are 



Generalized Deterministic Traffic 



Rules 3 




Figure 2: Fundamental diagram for lZ m .i for several different values of m. Vertical 
axis represents the flow <f) = pv. 
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moving with maximum velocity. In the "jammed phase" (above p c ), the average 
speed is equal to (1 — p)/p, just as for rule 184. The plot of the flow <fr = pv versus p 
for several different values of m, called the fundamental diagram, is shown in Figure 
It is clear that when the speed limit m increases, the total throughput of the 
highway increases. 



2. Monitored Traffic 

How can we increase the fluidity of the flow without increasing the speed limit? 
It is obvious that the source of inefficiency in rule 184 is the "defensive attitude" of 
the drivers. Consider, for example, the following configuration: 
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The first line represents locations of cars A, B, C at time t, and the second line 
their location at time t + 1. Zeros represent empty sites. Since driver A doesn't 
know whether car B is going to move or not, it is safer for him not to move. If he 
could see further than just one site forward, he could predict that car B will move 
forward, and that the site in front of him will be vacant at the next time step. 

A simple rule which incorporates such a "prediction" mechanism can be con- 
structed in a following way. Let us say that a driver at site i first checks whether in 
a block of k sites directly in front of him — i.e., from site i + 1 to site i + k — at least 
one site is empty (we will say that to be able to get this information cars have to 
be "monitored"). If it is, he moves his car by one site to the right, even if site i + 1 
is occupied since he knows that the car at site i + 1 will move because all drivers 
follow the same rule. If all sites from i + 1 to i + k are occupied, the car at site i 
does not move. Such a his will be denoted by Ri t k- For k = 2, the configuration 
discussed before will evolve as 
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As we see, cars A and B now move as one block, as long as the site in front of car B 
is empty. Such "blocking" significantly improves fluidity, as the fundamental dia- 
gram in Figure || demonstrates. We immediately notice that fundamental diagrams 
corresponding to rules 72-i,„ and 7Z n ,i are mirror images of each other with respect 
to the p = 0.5 line. Indeed Hi, n and lZ nt i are closely related: If cars are moving to 
the right according to rule lZi. n , then empty sites are moving to the left according 
to rule lZ n ,i. The critical density, which corresponds to a periodic configuration 
with a period consisting of k cars followed by one empty site equals p c — k/(k+ 1). 
Above p Cl a similar relationship holds, i.e., p c = kj(k + v). Hence 



1 if p< k/(k + l), 

k{l — p)/p otherwise. 



(4) 
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Figure 3: Fundamental diagram for 1Zi,k for several different values of k. 
3. Accidents 

Rules lZ n ,i and TZ\. n may be considered as two complementary driving strategies. 
Drivers obeying rule lZ n ,i drive with higher speed, but they keep the distance to the 
preceding car larger. Those who comply with rule TZi tU drive more slowly, however 
they keep a smaller distance between consecutive cars. Since the maximum possible 
flow for both rules is the same (it equals n/(n+ 1)), one could ask which rule is 
better in term of driving safety. Of course, if all drivers follow the same rule, either 
lZ n ,i or IZi^n, there is no problem , since no accident can occur. But, in the real 
world, some drivers are careless, have slower reaction time, defective brakes etc., 
and accidents do occur. On the single-lane highway, only one type of accident may 
occur: a car bumps into the preceding car which abruptly decreased its velocity. 
Thus, we can say that cars which decrease their velocity are potential causes of 
accidents. In our model, such potentially dangerous cars can be identified as cars 
which, at time t, have a smaller velocity that at time t—1 (by definition the velocity 
at time t is x(t) — x(t — 1), where x(t) is a position of a given car at time t). A 
simple mean-field estimate of a number of such cars can be carried out if we neglect 
time correlations between velocities, i.e., when we assume that v t -i and v t are not 
correlated. For rule 7£ TO ,i, the probability that a given car is slowing down will then 
be given by 

m i—1 

p(v t < «t_i) = J2 E p ^ = = o. (5) 

i=l j=0 
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where P(vt = j) denotes the probability that a car has velocity j at time t. Let 
us further assume that the length of the lattice is L, the number of cars is N, 
and the number of cars which have velocity i is iVj. In the stationary state JV< is 
time-invariant, thus we can write 

ra i—1 

p(v t <v t -i) = ^2Y1 n i n ^ ( 6 ) 

i=i j=o 

where n, = Ni/N. The values of Uj's for 7£ m ,i are given b 

n-k — n (l — n ) k if k < m (7) 
n m = (1 - n ) m 

Using these expressions, we obtain 

i—1 tn— 1 i—1 



P(u t < u t _i) =^nj-n OT + £ 

3=0 i=l 3=0 

m— 1 

= n (l -n ) m £(l-no) J + 



j'=o 

+ «o £ EC 1 -""W 1 -"»)*> 

i=l 3=0 

and after computing all sums, 

P(v t <v t - 1 ) = l-n -(l-n ) 2m + 

| (1 - n ) 2m + 2np - 1 - ng 
2 - rt 

As shown by Fukui and Ishibashil, when p > p c , no is a solution of the following 
nonlinear equation: 

(l-n )[l-(l-n ) m ] 1 



1, (8) 



n p 
If, for example, m = 2, this equation can be solved and uq is given by 



The fraction of cars which are slowing down equals 



3p 2 -5p-2 3 + p / 4~3^ 

p(ut < Vt - i] = — v — + v ~- (10) 

Below ,o c , P(«t < fj-i) is of course zero, since all cars are moving with constant 
speed m. 
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Figure 4: Fraction of cars which are slowing down for 72i,2 (□) and 72.2,1 (°). Solid 
and dashed lines represent the mean-field approximation for, respectively, rules 72i,2 
and 722, 1. 

For lZi t k, the mean-field approximation is much simpler, because all cars are 
either stopped or moving with speed 1 and v = X no + 1 X rc,i = m, thus 

P{v t < v t -i) = n ni = (1 - ni)m = (1 - v)v. (11) 

For k — 2, equation ([|) yields v = 2(1 — p)/p (above the critical density), therefore 

P(„ <„-,)= 1(1 - rt(: "- 2) . (12) 

P 

As before, below p c the fraction of cars which are slowing down is zero. 

The fraction of slowing cars obtained from computer simulations for both 72-2, 1 
and 72.1,2 is shown in Figure |], together with mean-field approximation curves given 
by equations ( |l0| ) and (|l2|). Although mean-field predictions overestimate P(vt < 
t>t_i), they are not very far from "experimental" results. One feature, however, is 
apparent: accident probability, which should be proportional to P(v t < ft-i), is 
much higher for 72.2,i than for 72i,2 if the density of cars is below 0.8. Above 0.8, rule 
72i,2 become more "dangerous", although the difference between rules diminishes 
as p approaches 1. It is also remarkable that in the case of rule 72i,2 we have no 
accidents up to p = 2/3, while rule 722, 1 becomes very dangerous quite fast, being 
worst at approximately p = 0.5. 
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4. Generalization 

Let us now construct a general rule which allows both higher speed limit m and 
blocking of A;-th order, with m and k both larger than 1. Consider a car located 
at site i. Its driver first locates the nearest gap (cluster of empty sites) in front of 
him whose length is g. If the first empty site (i.e., the first site belonging to the 
mentioned gap) is farther than i + fc, the car does not move, otherwise, it moves to 
site i + v site, where v — mm(g,m). We will refer to this rule as TZ m ^. We could 
now expect that, as before, the "perfect configuration" is a periodic sequence whose 
period consists of k cars followed by m empty sites. The critical density is then 
p c = k I (k + m), and this should hold for any p larger than p c , i.e., p = k/(k + v). 
Below p ci we expect v — m. For the average velocity, this leads to the following 
expression: 

v = 

Similarly, the flow is given by 

= 

Unfortunately, this simple reasoning has some flaws. Although in the vicinity of 
p = and p = 1 equations (|l^) and (|l4|) describe the behavior of lZ m ,k correctly, 
they fail to give a correct description in the intermediate region. This is clearly 
illustrated by the fundamental diagram, which, according to (|l4|), should be always 
tent-shaped, with a peak at p c = kj(k + m). In practice, however, the shape is 
quite different, as shown in Figure ^[ There is no single peak at p c , but instead, 
the flow stays at its maximum value over an extended interval of p values, and the 
fundamental diagram looks like "a tent with a flat roof" . This means that, unlike 
rules 1Zi t k an d Ti-m,i, rule lZ m .k for m, k > 1 does not exhibit a phase transition at 
p c . This is due to the fact that the "perfect configuration" with period consisting of 
k cars followed by m empty sites is no longer stable, like it was the case for rule IZi.k 
or IZm.i- In fact, the phase transition does occur, but it is, in a sense, spread over 
an interval. To be more precise, let us first note that v is a linear combination of 
n,-'s. Consequently, if, as functions of p, any of the rii has a discontinuous derivative, 
this will show up in v. Figure ^| shows an example of such a "velocity spectrum" , 
i.e., plots of rii(p) for i = 0,1,2,3. We can see that derivatives of all n^s have 
discontinuities at some point, but they occur at different values of p. The first 
transition occurs at p = 0.33 (n^ = 1 — > 723 ^ 1, Tlx = — > ni 7^ and ri2 = — > 
ri2 =/= ), whereas at p « 0.51 we have another transition uq — — > uq 7^ 0. For 
rule lZ m ,i, all these transitions occurred at the same point, and we consequently 
observed a single transition in v. 

Another interesting feature of the fundamental diagram of TZ m ,k is the maximum 
flow value. From (|l4| ) one would expect 4> m ax = km,/(k + m), which could be as 
large as we want if we choose the right m and k values. Figure f| however, clearly 
shows that this is not the case: there seems to be a cutoff at <fi slightly below 1. In 



m if p < Pc, , 13 ^ 

k{\ — p)/p otherwise. 



pm if p<p c , , u , 

(1 — p)k otherwise. 




fact, simulations performed with a wide array of k and m values suggest that the 
flow can never be larger that 1, regardless of k and m. For large k and m, the cutoff 
occurs almost exactly at 1, so the flow is well approximated by the piecewise linear 
function 

!pm if p < 1/m, 

1 if l/m< p < (k- l)/k, (15) 

(l-p)fe ifp > (fc- l)/fc. 

"Large" means equal or larger than 3, as even for m = k = 3 the above formula is 
fairly accurate (the cutoff occurs at 0.98, instead of 1). In any case, (f> = 1 is the 
absolute maximum flow in the stationary state. It is not possible to make it larger 
by increasing either k orm. 

5. Possible Modifications 

Driving strategies based on lZ m .k are, of course, not the only possible ways of 
increasing traffic fluidity. In the case of rule lZ m ,ki the decision of a driver where to 
move the car at the next time step is based only on the size of the nearest gap g nea r, 
located within the k + m — 1 sites in front of the car. Instead of considering the 
nearest gap, drivers could base their decisions on the largest gap. More precisely, 
let us say that g max is the largest continuous block of empty sites, entirely located 
between sites i and i + k + m. If, for example, m = 3, k = 7, consider the 
configuration . . . H)100I10000 . . ., the largest continuous block of zeros within m + 
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Figure 6: Velocity spectrum for 7^.3,2 . Spectrum computed after 5000 iterations, 
lattice size L = 2000. 

A; — 1 = 9 sites in front of the underlined car has a length g max — 3. We now define 
rule TZ m .k such that drivers move their cars by misx{g max , m) sites to the right. 
Since gmax > gnear, we can expect that the average speed (or flow) for 1Z m ,k should 
be larger or equal than the average speed (or flow) for TZ m ,k- To demonstrate that 
this is indeed the case, let us compare rules 7^.34 and Figure 0a shows that 

up to p w 0.43, 7^3,1 behaves exactly like 7^.34. When the density increases beyond 
this value, additional phase transitions occur. While rule 7^.3.1 exhibits just a single 
phase transition at p — 0.25, 1Zz t i exhibits an entire cascade, at p = 0.25, 0.43, 0.57, 
and 0.75. Interestingly, looking at the fundamental diagram (Figure 0b) we can see 
that it is fully symmetric with respect to p = 0.5. This means that the second peak 
has the same shape as the peak in fundamental diagram of (as a consequence 
of the duality between 72.3,1 and li-i.z). Therefore, 7^.34 can be considered as a sort 
of the "mix" of the aforementioned rules, 

v ~ J ^3,1 if p < 0.5, , 
n ^ ~ { H ll3 ifp>0.5. (16) 

The w sign reflects the fact that this relation is not true in the vicinity of 0.5. 

Another interesting detail which should be mentioned pertains to the maximum 
flow value. As in the lZ m .k case, we investigated the fundamental diagram of lZ m .k 
for many different m and k values, and it seems that also in this case 4> cannot 
exceed 1. Nevertheless, we found no simple explanation of this seemingly general 
feature. 
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Figure 7: (a) Plot of average velocity v as a function of density p for 7?-3 ; i (•) and 
7^3,1 (solid line), (b) Fundamental diagram of 0-3,1 (•) and TZs y i (solid line). 
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